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We explore the evolutionary dynamics of two games— the Prisoner's Dilemma and the Snowdrift Game- 
played within distinct networks (layers) of interdependent networks. In these networks imitation and 
interaction between individuals of opposite layers is established through interlinks. We explore an update 
rule in which revision of strategies is a biased imitation process: individuals imitate neighbors from the same 
layer with probability p, and neighbors from the second layer with complementary probability I — p. We 
demonstrate that a small decrease of p from p = I (which corresponds to forbidding strategy transfer 
between layers) is sufficient to promote cooperation in the Prisoner's Dilemma subpopulation. This, on the 
other hand, is detrimental for cooperation in the Snowdrift Game subpopulation. We provide results of 
extensive computer simulations for the case in which layers are modelled as regular random networks, and 
support this study with analytical results for coupled well-mixed populations. 

Network science has registered tremendous breakthroughs in the last years. Among these findings are the 
development of analytical and computational tools to study large real networks, as well as the develop- 
ment of more realistic models. These studies interpreted real networks mostly as single, isolated entities. 
However, the deeper understanding of complex networks is showing that in fact they are generally organized as 
networks of networks, and therefore the main focus of scientific research is shifting to multiplex and interdepend- 
ent networks^. 

Interdependent networks are organized in two or more layers (subnetworks), in which the functioning of a 
node in one layer depends on neighbor nodes in other layers. Multiplex networks^'^, on the other hand, are 
composed of a single type of node and several types ("colors") of links. One example of a real multiplex network is 
the transportation network of a city: nodes represent locations in the city, and links are of different type depending 
on the transportation system(s) that connect two given locations (metro, autobus, etc). Interdependent networks 
have already been applied in studies in diverse areas of science^. There is significant progress in the understanding 
of their percolation properties and robustness^'^"^^ Furthermore, interdependent networks were also applied in 
the evaluation of seismic risk^^'^^, in the emergence of creativity^^ and on the impact of such population structure 
on voting outcome^^. 

Network science^^"^^, together with evolutionary game theory^^ helps to understand how cooperative behavior 
can emerge and evolve in structured populations of selfish individuals. Cooperators contribute with a cost so that 
another individual can receive a certain benefit. Defectors on the other hand do not contribute, yet reap the benefit 
(which we assume to be larger than the cost). Interactions between individuals are traditionally modelled by one- 
shot symmetric two-person games, including the Prisoner's Dilemma (PD)^^'^^, the Snowdrift Game (SG)^^ or the 
Stag-Hunt game (SH)^^. In structured populations individuals only interact with their nearest neighbors, and it 
becomes possible for cooperators to escape exploitation by defectors by forming clusters in which they support 
each other. This process is known as network reciprocity^^'^^ . 

Population structure for games was traditionally modelled as a single, isolated network. Only very recently 
research has started to evaluate the impact of more complex, interdependent networks on the evolution of 
cooperation. To the best of our knowledge, all previous works have implemented the same game in all layers^^"^°. 
Different layers however can represent distinct environments, characterized by their game rules, which can 
interact or be coupled. Consider the example of two companies and the joint network of contacts between their 
employees. The code of conduct in each company may be different, depending on its internal organization and 
goals. Nevertheless, employees from one company may interact with employees from the other company, because 
they know each other personally or for instance when establishing a business transaction. By interacting with an 
acquaintance from the other company, each employee can acquire new working strategies, which can afterwards 



SCIENTIFIC REPORTS | 4 : 4436 | DOI: 1 0.1 038/srep04436 



1 





C 


D 


c 


RsG 


^SG 


D 




^SG 




c 


D 


C 


^PD 


^PD 


D 


TpD 


PpD 



SG layer: 




PD layer: 



! ! I 




SG 

quadrant 



PD 

quadrant 



1 Tpo = Tse 2 



Figure 1 | Scheme of the population structure. The population is organized in two layers. Individuals establish intralinks with neighbors of the same layer, 
and interlinks with neighbors of the opposite layer. Depending on the layer in which they are located, individuals compute their payoff taking into 
account the PD payoff matrix or the SG payoff matrix, indicated on the left of the figure. On the right, we schematize the relation between the parameters 
of the two games, which we use in this paper. We assume that Tpo = Tsg = T, Ssg = S and Spo = — 5, with S G [0, 1], T G [1,2]. 



be imitated by their contacts inside each company. Another natural 
example is interaction between subpopulations of different cultural 
backgrounds. 

Here we investigate how imitation and interaction between indi- 
viduals playing different games influence the final cooperation levels. 
For that, we consider a population organized in two layers of equal 
size, in which one distinct game is played in each of the layers, as 
schematized in Fig. 1 (see also Methods for details). Nodes establish 
intralinks with neighbors of the same layer, and interlinks with neigh- 
bors of the opposite layer. In the following, we start by studying by 
computer simulations the case in which both layers are regular ran- 
dom graphs. Afterwards, we also consider the limit case of a well- 
mixed population of infinite size, also organized in two layers. We 
introduce bias in imitation: individuals imitate a neighbor from the 
same layer with probability p, or a neighbor of the opposite layer with 
probability (1 — p). We investigate the impact of varying p on the 
final level of cooperation reached in each of the layers (see Methods). 

Results 

Regular random layers. For reference, we present in Fig. 2 the final 
fraction of cooperators for the PD and the SG played in isolated 
populations. Panel A) shows the results of computer simulations in 
regular random networks, while panel B) shows the stationary 
solutions of the replicator equation x = x{l—x)\fc{x)—fo{x)] for 
an infinite well-mixed population, that is a population in which 
every individual has the same probability of interacting with 
anyone else^\ In the replicator equation, x represents the fraction 
of cooperators in the population, fcM and /d(^) stand for the 
mean fitness of cooperators and defectors respectively (with fdx) 
= xR + {I - x)S and /d = + (1 - x)P). The PD is the 
harshest social dilemma for cooperation, and therefore only a small 
region characterized by S close to 0 and T close to 1 is able to escape 
full defection, in regular random networks. The SG, on the other 
hand, is a coexistence game, in which the final fraction of 
cooperators x in an infinite well-mixed population is given by x = 
{P - S)/{R - S - T + P). 

To evaluate the impact of an interdependence between the PD and 
the SG on the evolution of cooperation, we start by considering a 
population in which each of the layers is modeled by a regular ran- 
dom graph with nodes of degree k = 4. This graph has no finite loops 
in the infinite size limit. We assume that each node establishes one 
interlink with a uniformly randomly chosen node from the opposite 
layer. That is, individuals engage in four interactions with neighbors 
of the same layer and one interaction with a neighbor from the other 
layer. 

The contour plots in Fig. 3 show the result of our simulations, 
namely, the final fraction of cooperators in each of the layers as a 



function of the game parameters T and S, for several biased imitation 
probabilities p (0 < p < 1). Whenp = 1 individuals can only imitate 
neighbors from the same layer, and the results obtained are qualita- 
tively similar to those known for isolated populations (see panel A) of 
Fig. 2). 

When p < 1 strategy transfer between layers becomes possible. 
Comparing the results for the PD layer forp = 0.9, p = 0.95 andp = 
1, we reach an interesting conclusion: although the fraction of coop- 
erators is still low, even a small probability of imitating a neighbor 
playing a distinct social dilemma (in this case, the SG) is sufficient to 
prevent full defection in the majority of the PD quadrant. 

Asp further decreases, for a given pair of game parameters T, S the 
fraction of cooperators in the stationary state becomes approxi- 
mately equal in both layers (note that Ssg — Spo = —S with S 
e [0, 1]). However, a careful inspection of the contour plots pre- 
sented in Fig. 3 for 0.1 < p < 0.5 shows that the fraction of coopera- 
tors in both layers is not exactly identical (i.e. the contours are not 
exactly symmetrical): for a fixed pair of T and S values, the PD layer 
tends to have a slightly lower fraction of cooperators. 

Finally, the results for p = 0 are surprisingly different from the 
ones for higher p. This contrast indicates the different structures of 
absorbing states at different values of p. Whenp = 1 four absorbing 
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Figure 2 | Final fraction of cooperators in single, isolated populations. 

The contour plot of panel A) shows the final average fraction of 
cooperators for the PD and SG when these social dilemmas are played in an 
isolated regular random network of size N = 10^ and node degree k = A, 
Panel B) shows the corresponding results for an infinite well-mixed 
population, solution of the replicator equation x = x{l—x)\fc {x) — /d (x)] , 
where ^(x) and^(x) stand for the mean fitness of cooperators and 
defectors respectively^^ Red corresponds to full defection, blue 
corresponds to full cooperation. The parameter P is set to 1.0 for panel A). 
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Figure 3 | Final fraction of cooperators when layers are modelled as random regular graphs, for various values of probability p. Contour plots show the 
final average fraction of cooperators on both the PD (S < 0) and SG (S > 0) layers as a function of the game parameters Tand S, starting from 50% of 
cooperators and defectors randomly distributed in each of the layers. The dashed line separates the results for the SG layer (upper quadrants) from 
the results for the PD layer (lower quadrants). Red corresponds to full defection, blue corresponds to fijU cooperation (x stands for the fraction of 
cooperators in a given layer). 



states exist, and both layers reach homogeneous states: both reaching 
full defection (full cooperation), or one layer reaching full coopera- 
tion and the other full defection (and vice versa). When 0 < p < 1, 
individuals can imitate neighbors from both layers, and therefore 
only the two absorbing states in which both layers reach the same 
homogeneous state are possible. In contrast, when p = 0, each indi- 
vidual can only imitate his neighbor connected through the interlink. 
As a result, there is an infinite number of absorbing states. When all 
interlinks connect pairs of individuals with equal strategies, evolu- 
tion stops, as schematized in panel A) of Fig. 4. Since strategies and 
interlinks are distributed randomly in the population, on average less 
than N/2 interlinks will actually connect individuals of different 
strategies, and the population will rapidly evolve to an absorbing 
state. 

Figure 5 shows the cross section of the contour plots of Fig. 3 at 
fixed S = 0.5 for several p values. This illustrates the different beha- 
vior of the level of cooperation in each of the layers for decreasing/). 
Excluding the results for p = 0 (which are a direct consequence of 
allowing exactly one interlink for each node) we observe that the 
fraction of cooperators monotonously decreases for decreasing p in 
the SG layer, while in the PD layer it monotonously increases. As we 
will see in the following, this behavior will not be maintained when 
this population structure is replaced by a well-mixed population. 
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Figure 4 | Absorbing states for p= 0. Each panel represents a snapshot of 
the two layers seen in profile. Panel A) schematizes one possible absorbing 
state when layers are modelled as regular random networks. Since each 
individual can only imitate his neighbor connected through the interlink, 
the population reaches an absorbing state organized as an arbitrary set of 
independent pairs (cooperator-cooperator and defector-defector) in 
which one member is in one layer and the other is in the other layer. Panel 
B) represents one of the two possible absorbing states in a well-mixed 
population. Since each node establishes interUnks to every node from the 
opposite layer, evolution ceases when both layers reach full cooperation or 
full defection. 



Well-mixed population. We now proceed to analyze the general case 
of a population modelled by a fully connected network organized in 

two layers. Each layer has N/2 nodes and -{N/2){N/2 — l) 

intraUnks. Contrary to the previous case of regular random layers, 
in which each node established one interlink, we now assume that 
each node establishes N/2 interUnks to the N/2 nodes of the opposite 
layer, in accordance with the structure of a fully connected network. 
Therefore, each node has N — 1 neighbors: N/2 — 1 in the same layer 
as itself, and N/2 in the opposite layer. Up^ and Usg stand for the 
number of cooperators in the PD (SG) layer, respectively. In this case, 
mean-field theory is exact and the evolution of population can be 
expressed completely by (t) and Hsg (t)- 

We obtain the solutions for the stationary state of this model in the 
infinite population limit, N^ ^. For that, we adopted the variables 
^PD = npY)/{N/2) and x^g = /1sg/(^/2) (xpd ^ [0, 1] and Xsg ^ 
[0, 1]). The mean fitness of individuals in each of the layers is 
expressed in terms of the game parameters and the fraction of coop- 
erators and defectors. In the PD layer, the fitnesses of cooperators 




T 



Figure 5 | Final fraction of cooperators versus Tat fixed 5 = 0.5, for both 
layers. Results correspond to those presented in Fig. 3 when S = 0.5. The 
fraction of cooperators for the SG layer is shown in the upper plot, while 
the corresponding curves for the PD layer are presented in the bottom plot. 
For simplicity, we include only some of the p values studied in Fig. 3. 
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Figure 6 | Final fraction of cooperators on each layer of the well-mixed population, for various values of probability p. Contour plots show the 
final fraction of cooperators on both the PD (S < 0) andSG (S> 0) layers (solution ofEqs. (5) and (6)) as afijnction of the game parameters Tand S. The 
dashed line separates the results for the SG layer (upper quadrants) from the results for the PD layer (lower quadrants). Red corresponds to hall 
defection, blue corresponds to full cooperation (x stands for the fraction of cooperators in a given layer). Parameters: P = 1.0. 



(/c,pd) and defectors (/d,pd) are expressed as: 

/C,PD =^PD - (1 -Xpd)S + xsg - (1 -xsg)S, (1) 

/d,pd=^pdT + XsgT. (2) 

Similarly, in the SG layer the fitnesses of cooperators (/c,sg) and 
defectors (/d,sg) are 

/c,SG =^SG + (1 -^sg)S + Xpd + (1 -Xpd)S, (3) 

/d,sg = ^sg^ + ^pdT. (4) 



We obtain the following rate equations for the evolution of the 
fraction of cooperators in the PD and SG layers: 



:^PD = (1 -:^pd) \pXpDF(fc,PD -/d,pd) + (1 -p)xsGF(fc,SG -/d,pd)] 

-XpD [p(l - a:pd)^'(/d,pd -/c,pd) + (1 -p) (1 -XsG)F(fD,SG -few)]. 



(5) 



^SG = (1 -^Sg) \pXsGF(fc,SG -/d,Sg) + (l -p)^PDi^(/c,PD -/d,Sg)] 

-Xsg[p(1 -XsG)F(fD,SG -/cSg) + (l -p)(l -XpB)F(fD,FD -/c,Sg)] 



(6) 



For instance, in Eq. (5) the factor (1 — Xpd) stands for the prob- 
ability of selecting a defector from the PD layer. The second factor, in 
the square brackets, accounts for the two possible scenarios: with 
probability p the defector can imitate a cooperator from the PD layer, 
while with complementary probability (I — p) the defector can imit- 
ate a cooperator that belongs to the SG layer. The interpretation is 
analogous for the second line of this equation, as well as for Eq. (6). 
The stationary state solution corresponds to that obtained from Eqs. 
(5) and (6) when = 0 and i:sG = 0. 

Figure 6 shows the final stationary states of this mean-field model 
for different values of probability provided by the numerical solu- 
tions of Eqs. (5) and (6) for p = 1.0. These plots actually explain our 
simulations in the previous section for coupled regular random net- 
works (see Fig. 3). One can clearly see that Figs. 3 and 6 show qual- 
itative similarities, with only one dramatic exception, p = 0, which 
will be discussed below. In particular, in the case ofp= 1 both Figs. 3 
and 6, as well as Figs. 5 and 7, are close to the corresponding results 
for isolated populations (random regular and well-mixed, respect- 
ively, compare with Fig. 2). 

One can see from Fig. 6 that a small probability of imitating a 
neighbor from the SG layer is beneficial for the levels of cooperation 



reached in the PD layer. The fmal fraction of cooperators in each of 
the layers becomes approximately equal as p decreases, yielding 
almost symmetrical contours for 0.1 < p < 0.5. However, the frac- 
tion of cooperators is in general slightly lower in the PD layer than in 
the SG layer for a fixed pair of game parameters T and S. 

Note some difference between Figs. 3 and 6 at sufficiently low, still 
non-zero, p. With diminishing p, the region of full defection in 
coupled regular random populations (Fig. 3) approaches some limit. 
In contrast, for decreasing/), coupled well-mixed populations (Fig. 6) 
do not show any area with full defection, but the region with high 
level of defection becomes wider and wider asp decreases. Finally, at 
p = 0, the full defection state is realized in the entire T — S plane. This 
difference is also observed between Figs. 5 and 7. Especially it man- 
ifests itself in two different behaviors of fractions of cooperators in 
the PD layer in Figs. 5 and 7, monotonous with p and non-mono- 
tonous, respectively. 

At p = 0 the difference between the coupled regular random net- 
works and well-mixed population is dramatic. While for the coupled 
well-mixed populations, the final state is full defection for the whole 
T — S plane (Fig. 6), for regular random networks we observed 




p = 0.9 
p = 0.5 
p = 0.25 
p = 0.1 



Figure 7 | Final fraction of cooperators versus Tat fixed 5 = 0.5 for both 
layers of a fully connected network. Results correspond to those presented 
in Fig. 6 when S = 0.5. The fraction of cooperators for the SG layer is shown 
in the upper plot, while the corresponding curves for the PD layer are 
shown in the bottom plot. 
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coexistence of defection and cooperation for the whole T — S plane. 
The reason for this controversy is the fundamental difference 
between the structures of the absorbing states in these two systems. 
Clearly, in the mean-field system (well-mixed) there are two absorb- 
ing states: full cooperation and full defection. For p = 0, the full 
defection scenario is realized. In contrast, in the coupled regular 
random populations, at p = 0 the system has an infinite number of 
absorbing states. We already explained that this state is organized as 
an arbitrary set of independent pairs (cooperator-cooperator and 
defector- defector) in which one member is in one layer and the other 
is in the other layer. The final fraction of cooperators and defectors in 
the population for a given point on the T — S plane is determined by a 
specific initial condition. 

It is interesting to relate our findings with observations in other 
works^^''*^ which studied the role of biased fitness (utility) on the level 
of cooperation in interdependent structured populations. Although 
these works explored a bias and network architectures very different 
from our study, they reported the existence of an optimal control 
parameter value for cooperation, similarly to what we notice for a 
well-mixed population. 

We can also relate our model to the multi-network model^^'^^, 
where links between nodes are of two sorts: the interaction links 
through which individuals play and the imitation links through 
which individuals revise their strategies. All links in our model are 
used both for interaction and imitation. The principal difference is 
that interaction occurs with probability 1 for every link, while imita- 
tion takes place with probability p along intralinks and the comple- 
mentary probability I — p along interlinks. 

Discussion 

We explored the evolution of cooperation in two games— the 
Prisoner's Dilemma and the Snowdrift Game — played within an 
interdependent network. We demonstrated that as the probability 
to imitate neighbors from the opposite layer increases, the final level 
of cooperation in each of the layers behaves differently. A small 
probability 1 — p is always detrimental for the level of cooperation 
in the SG layer, regardless of the population structure (well-mixed or 
when layers are modelled as regular random networks). In the case of 
a well-mixed population we observed an intermediate optimal value 
ofp for which the fraction of cooperators in the steady state is max- 
imal in the PD quadrant. On the other hand, when layers are regular 
random graphs, the level of cooperation in the PD layer registers a 
significant increase from full defection as p decreases from 1, as is 
shown in Fig. 3. We have also shown that different arrangements of 
interconnections between layers result in fundamentally different 
structures of absorbing states forp = 0, that is the limit case in which 
individuals can only imitate neighbors through the interlinks. 

The constraint between the game parameters of the PD and the SG 
implemented in this work is one of possible options. Based on our 
mean-field approach we have verified that our conclusions remain 
qualitatively valid for other relationships between game para- 
meters^^. In this work we studied a representative case of P = 1. In 
the Supplementary Information we consider other values of j]. 

The Prisoner's Dilemma and the Snowdrift Game are two of the 
possible two-person games. Future research directions include estab- 
lishing an interdependence between other two- and AT-person games, 
as well as extending the biased imitation process here reported to the 
case in which layers are strongly heterogeneous, such as scale-free 
networks^^. 
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in which the entries represent the payoff earned by the row player. When both 
individuals opt for cooperation (defection), they both receive a Reward, R 
(Punishment, P). When a cooperator meets a defector, the cooperative player received 
the Sucker's payoff (S), while the defective player receives the Temptation (T) for 
defecting. 

The relative ordering of the game parameters R, S, T and P defines four different 
two-person games^^. In this study, we consider the Prisoner's Dilemma (PD) (char- 
acterized by r > > P > Sy^ and the Snowdrift Game (SG) (characterized by T > P 

> S > Py\ We fix P = 1, P = 0 as is conventional"'. 

Population structure. We consider a population organized in two layers 
(schematized in Fig. 1), in which each node represents a different individual. In each 
layer all individuals play the same game, the PD in one layer and the SG in the other. 
We refer to links between nodes of opposite layers as interlinks, and links between 
nodes of the same layer as intralinks. In the absence of interlinks, each layer stands as 
an isolated population in which one single game is played. When layers are modelled 
as regular random graphs, we assume that each node establishes one interlink to a 
randomly chosen node from the opposite layer. On the other hand, when studying a 
two-layer fully connected network, we assume that each node establishes interlinks to 
each of the nodes in the opposite layer. 

Relation between game parameters. We assume that individuals in the PD layer 
adopt the game parameters Ppo, Spo, Tp^ and PpD (with Tp^ > Rp^ > Ppo > Spo). 
Similarly, individuals in the SG layer consider Psg> 5sg> ^sg and Psg (with Tsg > -^sg 

> SsG > ^'sg)- The game parameters for the PD and the SG can be related in several 
ways. We opt for Rpi^ = Psg = 1> Ppd = -Psg = 0 and Tpo = Tsg = Twith T ^ [1,2]. 
Furthermore, we constrained Ssg = S and SpD = —S, with S e [0, 1]. The resulting 
payoff matrix for the SG is 



C D 

C / 1 s 
D\T 0 

while the payoff matrix for the PD layer is given by 

C D 
C/l -S^ 
D\T 0 



(8) 



(9) 



The difference between matrix (8) for the SG layer and matrix (9) for the PD layer is in 
the payoff obtained by a cooperator when interacting with a defector, the sucker's 
payoff. When interacting with a defector, a cooperator located in the SG layer earns S, 
while a cooperator located in the PD layer earns — S. We stress that this is only one of 
possible relations between the parameters of the two games. 

Individuals play the same game with all neighbors, including the neighbor(s) from 
the other layer. An individual's fitness represents a measure of his success, and 
corresponds to the payoff accumulated over all interactions in which he engages in 
one time-step. Through the interlinks, individuals obtain their payoff by taking into 
account the strategy of their co-player and the payoff matrix used in their own layer. 

Update dynamics. Individuals revise their strategies through a biased imitation 
process: with probability p they imitate a neighbor from the same layer, while with 
probability (1 — p) individuals imitate a neighbor from the opposite layer. Update is 
asynchronous: at each time-step we randomly choose one of the layers, and from it 
one randomly chosen individual. A, to revise his strategy. Then we choose with 
probability p a neighbor from the same layer, and with complimentary probability 
(1 — p), a neighbor from the opposite layer. If A and his chosen neighbor, B, have 
different strategies, A imitates the strategy of B with a probability F (/b — /a) 
proportional to their fitness difference and given by the Fermi distribution*^'^*^. 



1 



l+e-W.-/*) 



(10) 



/a (/b) stands for the fitness of A (B), and represents the intensity of selection, which 
regulates the accuracy of the imitation process. In the limit ^ 0 evolution proceeds 
by random drift, while fi^^ corresponds to the deterministic limit of pure copying 
dynamics: A imitates B iff/e >/a, with probability 1. We consider the respresentative 
case of B = 1.0. 



Methods 

Evolutionary games. We implement standard two-person games. Each individual 
can act as an unconditional cooperator (C) or as an unconditional defector (D). The 
possible payoffs resulting from the interaction between two players are resumed in the 
payoff matrix 



Simulations. The results of Figs. 3 and 5 were obtained for a population of total size N 
= 2000, in which both layers had an equal size, N/2 = 1000. Each point in the contour 
plots corresponds to an average over 10^ trials. In each trial, we start by randomly 
selecting one over 10^ realizations of a regular random network to model each layer, 
and establish N/2 interlinks between the layers. Each interlink is established by 
uniformly randomly choosing two end nodes, one from each layer, allowing each 



SCIENTIFIC REPORTS | 4 : 4436 | DOI: 1 0.1 038/srep04436 



5 



node to establish exactly one interlink. Each trial starts by randomly distributing 50% 
of cooperators and defectors in each of the layers. We allow the population to evolve 
for lO'^ generations (where 1 generation corresponds to Nil strategy revisions), after 
which we average the fraction of cooperators over the next 10^ generations. 



1. Buldyrev, S. V., Parshani, R., Paul, G., Stanley, H. E. & Havlin, S. Catastrophic 
cascade of failures in interdependent networks. Nature 464, 1025-1028 (2010). 

2. Nicosia, V., Bianconi, G., Latora, V. & Barthelemy, M. Growing multiplex 
networks. Phys. Rev. Lett. Ill, 058701 (2013). 

3. Cozzo, E., Banos, R. A., Meloni, S. & Moreno, Y. Contact-based social contagion in 
multiplex networks. Phys. Rev. E 88, 050801 (2013). 

4. Helbing, D. Globally networked risks and how to respond. Nature 497, 51-59 
(2013). 

5. Baxter, G. J., Dorogovtsev, S. N., Goltsev, A. V. & Mendes, J. F. F. Avalanche 
collapse of interdependent networks. Phys. Rev. Lett. 109, 248701 (2012). 

6. Parshani, R., Buldyrev, S. V. & Havlin, S. Interdependent networks: reducing the 
coupling strength leads to a change from a first to second order percolation 
transition. Phys. Rev. Lett. 105, 048701 (2010). 

7. Son, S.-W., Bizhani, G., Christensen, C, Grassberger, P. & Paczuski, M. 
Percolation theory on interdependent networks based on epidemic spreading. 
EPL 97, 16006 (2012). 

8. Gao, J., Buldyrev, S. V., Havlin, S. & Stanley, H. E. Robustness of a network of 
networks. Phys. Rev. Lett. 107, 195701 (2011). 

9. Zhou, D., Stanley, H. E., D'Agostino, G. & Scala, A. Assortativity decreases the 
robustness of interdependent networks. Phys. Rev. E 86, 066103 (2012). 

10. Li, W., Bashan, A., Buldyrev, S. V., Stanley, H. E. & Havlin, S. Cascading failures in 
interdependent lattice networks: The critical role of the length of dependency 
links. Phys. Rev. Lett. 108, 228702 (2012). 

11. Huang, X., Gao, J., Buldyrev, S. V., Havlin, S. & Stanley, H. E. Robustness of 
interdependent networks under targeted attack. Phys. Rev. E 83, 065101 (2011). 

12. Vespignani, A. Complex networks: The fragility of interdependency. Nature 464, 
984-985 (2010). 

13. Radicchi, F. & Arenas, A. Abrupt transition in the structural formation of 
interconnected networks. Nat. Phys. 9, 717-720 (2013). 

14. Poljansek, K., Bono, F. & Gutierrez, E. Seismic risk assessment of interdependent 
critical infrastructure systems: The case of european gas and electricity networks. 
Earthquake Eng. Struct. Dynam. 41, 61-79 (2012). 

15. Duenas-Osorio, L., Craig, J. I. & Goodno, B. J. Seismic response of critical 
interdependent networks. Earthquake Eng. Struct. Dynam. 36, 285-306 (2007). 

16. Csermely, P. The appearance and promotion of creativity at various levels of 
interdependent networks. Talent Development Excellence 5, 115-123 (2013). 

17. Halu, A., Zhao, K., Baronchelli, A. & Bianconi, G. Connect and win: The role of 
social networks in political elections. EPL 102, 16002 (2013). 

18. Dorogovtsev, S. N. & Mendes, J. F. Evolution of networks. Adv. Phys. 51, 
1079-1187 (2002). 

19. Dorogovtsev, S. N. & Mendes, J. F. Evolution of Networks: From Biological Nets to 
the Internet and WWW (Oxford University Press, 2003). 

20. Albert, R. & Barabasi, A.-L. Statistical mechanics of complex networks. Rev. Mod. 
Phys. 74, 47 (2002). 

21. Sigmund, K. The Calculus of Selfishness (Princeton University Press, 2010). 

22. Axelrod, R. & Hamilton, W. D. The evolution of cooperation. Science 211, 
1390-1396 (1981). 

23. Rapoport, A. Prisoners Dilemma: A Study in Conflict and Cooperation (University 
of Michigan Press, 1965). 

24. Sugden, R. The Economics of Rights, Co-operation and Welfare (Basil Blackwell 
Oxford, 1986). 

25. Skyrms, B. The Stag Hunt and the Evolution of Social Structure (Cambridge 
University Press, 2004). 

26. Nowak, M. A. & May, R. M. Evolutionary games and spatial chaos. Nature 359, 
826-829 (1992). 

27. Nowak, M. A. Five rules for the evolution of cooperation. Science 314, 1560-1563 
(2006). 

28. Santos, F. C. & Pacheco, J. M. Scale-free networks provide a unifying framework 
for the emergence of cooperation. Phys. Rev. Lett. 95, 098104 (2005). 

29. Santos, F. C, Santos, M. D. & Pacheco, J. M. Social diversity promotes the 
emergence of cooperation in public goods games. Nature 454, 213-216 (2008). 



30. Szabo, G. & Path, G. Evolutionary games on graphs. Phys. Rep. 446, 97-216 
(2007). 

31. Hauert, C. & Doebeli, M. Spatial structure often inhibits the evolution of 
cooperation in the snowdrift game. Nature 428, 643-646 (2004). 

32. Ohtsuki, H., Hauert, C, Lieberman, E. & Nowak, M. A. A simple rule for the 
evolution of cooperation on graphs and social networks. Nature 441, 502-505 

(2006) . 

33. Roca, C. P., Cuesta, J. A. & Sanchez, A. Effect of spatial structure on the evolution 
of cooperation. Phys. Rev. E 80, 046106 (2009). 

34. Gomez- Garderies, J., Reinares, I., Arenas, A. & Fiona, L. M. Evolution of 
cooperation in multiplex networks. Sci. Rep. 2, 620 (2012). 

35. Wang, Z., Szolnoki, A. & Perc, M. Interdependent network reciprocity in 
evolutionary games. Sci. Rep. 3, 1183 (2013). 

36. Gomez- Gar denes, J., Gracia-Lazaro, C, Floria, L. M. & Moreno, Y. Evolutionary 
dynamics on interdependent populations. Phys. Rev. E 86, 056113 (2012). 

37. Wang, B., Chen, X. & Wang, L. Probabilistic interconnection between 
interdependent networks promotes cooperation in the public goods game. /. Stat. 
Mech. 2012, PI 1017 (2012). 

38. Jiang, L.-L. & Perc, M. Spreading of cooperative behaviour across interdependent 
groups. Sci. Rep. 3, 2483 (2013). 

39. Wang, Z., Szolnoki, A. & Perc, M. Optimal interdependence between networks for 
the evolution of cooperation. Sci. Rep. 3, 2470 (2013). 

40. Wang, Z., Szolnoki, A. & Perc, M. Evolution of public cooperation on 
interdependent networks: The impact of biased utility functions. EPL 97, 48001 
(2012). 

41. Taylor, P. D. & Jonker, L. B. Evolutionary stable strategies and game dynamics. 
Math. Biosci. 40, 145-156 (1978). 

42. Ohtsuki, H., Nowak, M. A. & Pacheco, J. M. Breaking the symmetry between 
interaction and replacement in evolutionary dynamics on graphs. Phys. Rev. Lett. 
98, 108106 (2007). 

43. Ohtsuki, H., Pacheco, J. M. & Nowak, M. A. Evolutionary graph theory: breaking 
the symmetry between interaction and replacement. /. Theor. Biol. 246, 681-694 

(2007) . 

44. Santos, M. D., Dorogovtsev, S. N. & Mendes, J. F. F. in preparation. 

45. Barabasi, A.-L. & Albert, R. Emergence of scaling in random networks. Science 
286, 509-512 (1999). 

46. KoUock, P. Social dilemmas: The anatomy of cooperation. Annu. Rev. Sociol. 
183-214 (1998). 

47. Traulsen, A., Nowak, M. A. & Pacheco, J. M. Stochastic payoff evaluation increases 
the temperature of selection. /. Theor. Biol. 244, 349-356 (2007). 

48. Traulsen, A., Pacheco, J. M. & Nowak, M. A. Pairwise comparison and selection 
temperature in evolutionary game dynamics. /. Theor. Biol. 246, 522-529 (2007). 



Acknowledgments 

This work was partially supported by the FCT project PTDC/MAT/1145 15/2009 and the 
FET proactive IP project MULTIPLEX number 317532 and the fellowship grant SFRH/ 
BPD/90936/2012. 



Author contributions 

M.D.S., S.N.D. and J.F.F.M. equally designed and did the research, discussed the results and 
wrote the paper together. 

Additional information 

Supplementary information accompanies this paper at http://www.nature.com/ 
scientificreports 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Santos, M.D., Dorogovtsev, S.N. & Mendes, I.F.F. Biased imitation 
in coupled evolutionary games in interdependent networks. Sci. Rep. 4, 4436; DDI: 10. 1038/ 
srep04436 (2014). 



©0®© 



This work is licensed under a Creative Commons Attribution- 
NonCommercial-NoDerivs 3.0 Unported license. To view a copy of this license, 



visit http://creativecommons.Org/licenses/by-nc-nd/3.0 



SCIENTIFIC REPORTS | 4 : 4436 | DOI: 1 0.1 038/srep04436 



6 



